plot(t,f(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal f(t)=t^2')
stem(n,g(n),'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal f(n)=n^2 ')
plot(t,f(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal f(t)=e^{-t} ')
y = piecewise(t < 0,0,t >= 0,exp(-t))
y =

fplot(y,[-3 5], 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal causal f(t)')
t = -5:0.01:3; % ajuste de t
plot(t,f(t+2), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal f(t+2)')
plot(t,f(t)+2, 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal f(t)+2')
t = 0:0.01:8; % ajuste de t
plot(t,f(t-3), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal f(t-3)')
plot(t,f(t)-1, 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal f(t)-1')
plot(t,f(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-5 5 -0.2 1.2])% ajuste visualización
title('Gráfica de la señal f(t)')
t = -pi:0.01:0; % ajuste t
plot(t,f(-t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-5 5 -0.2 1.2])% ajuste visualización
title('Gráfica de la señal f(-t)')
t = 0:0.01:pi/2; % ajuste t
plot(t,f(2*t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-5 5 -0.2 1.2])% ajuste visualización
title('Gráfica de la señal f(2t)')
t = -pi/2:0.01:0; % ajuste t
plot(t,f(-2*t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-5 5 -0.2 1.2])% ajuste visualización
title('Gráfica de la señal f(-2t)')
t = 0:0.01:3*pi; % ajuste t
plot(t,f((1/3)*t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-10 10 -0.2 1.2])% ajuste visualización
title('Gráfica de la señal f((1/3)t)')
t = -3*pi:0.01:0; % ajuste t
plot(t,f((-1/3)*t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-10 10 -0.2 1.2]) % ajuste visualización
title('Gráfica de la señal f((-1/3)*t)')
plot(t,f(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-5 5 -2.2 2.2])% ajuste visualización
title('Gráfica de la señal f(t)')
plot(t,-2*f(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-5 5 -2.2 2.2])% ajuste visualización
title('Gráfica de la señal -2f(t)')
plot(t,2*f(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-5 5 -2.2 2.2])% ajuste visualización
title('Gráfica de la señal 2f(t)')
plot(t,(-1/3)*f(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-5 5 -2.2 2.2])% ajuste visualización
title('Gráfica de la señal (-1/3)f(t)')
plot(t,(1/3)*f(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-5 5 -2.2 2.2])% ajuste visualización
title('Gráfica de la señal (1/3)f(t)')
plot(t,-f(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
axis([-5 5 -2.2 2.2])% ajuste visualización
title('Gráfica de la señal -f(t)')
g(t) = abs(t)*exp(-abs(t))
g(t) = 
g(-t)
ans = 
h(t)= t^3*cos(t)
h(t) = 
h(-t)
ans = 
f(t) = cos(t)+sin(3*t)-1
f(t) = 
f(-t)
ans = 
disp('f no es impar ni par')
f_p(t) = (f(t) + f(-t))/2
f_p(t) = 
f_p(-t)
ans = 
f_i(t) = (f(t) - f(-t))/2
f_i(t) = 
f_i(-t)
ans = 
fplot(f_p(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal f_p(t)')
fplot(f_i(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal f_i(t)')
fplot(f(t), 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal f(t)')
f(t) = exp(-t)
f(t) = 
g(t) = f(t)*u(t)
g(t) = 
fplot(g(t),[-3 5], 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal causal g(t)')
f(t) = u(t+1) - u(t-1) - 2 * u(t);
fplot(f(t),[-1.5 1.5], 'LineWidth',2)
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Gráfica de la señal causal g(t)')
plot3(t,real(f),imag(f),'LineWidth',2)
title('Gráfica de la señal f(t)=e^{st} ')
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('Parte imaginaria y parte real')
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('tiempo y parte real')
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
title('tiempo y parte imaginaria')
syms s t y(t) x(t) Y(s) X(s) yy
Li= diff(y(t),t,2)+ 5* diff(y(t),t) + 6*y(t)
Li =

Ld = diff(x(t),t)+x(t)
Ld =

ecu = Li == Ld
ecu =

LI = laplace(Li)
LI =

LI = subs(LI,[laplace(y(t), t, s), y(0),subs(diff(y(t), t), t, 0)],[Y(s),2,1])
LI = 
LD = laplace(Ld)
LD = 
LD = subs(LD,[laplace(x(t), t, s),x(0)],[X(s),0])
LD = 
ecu = LI==LD
ecu = 
ecu = subs(ecu,X(s),laplace(exp(-4*t)*heaviside(t)))
ecu =

Y(s)=simplify(solve(ecu,yy))
Y(s) =

Y(s)= partfrac(Y(s))
Y(s) =

y(t) = ilaplace(Y(s))*heaviside(t)
y(t) =

function laplace2016a(a,b,ciy,xi,t0)
% a coeficientes de las derivadas de la salida menor a mayor [a_0, ..., a_n]
% b coeficientes de las derivadas de la entrada menor a mayor [b_0, ..., b_m]
% ciy condiciones iniciales de la salida de menor a mayor [y(0), y(0)^(n-1)]
% xi función de entrada en terminos de la variable simbolica t previamente
% declarada en el command window
% t0 tiempo final para graficar la solucion, la derivada, y la segunda
% ejemplo: resolver y^(3)+y^(2)+2y^(1)+2y=3x^(2)-x^(1)+2x con y^(2)(0)=1 y^(1)=3
% y(0)=2, x(t)=exp(-t)cos(t)u(t), para 10 segundos, se resuleve como
% laplace2016a([2 2 1 1],[2 -1 3],[2 3 1],exp(-t)*cos(t)*heaviside(t),10)
syms y(t) Y(s) x(t) X(s) Yy fp;
edd=edd+a(i)*s^(i-1)*Y(s);
edd=edd-a(i)*(s^(i-1-k)*ciy(k));
edi=edi+b(i)*s^(i-1)*X(s);
% edi=edi-b(i)*(s^(i-1-k)*cix(k));
mensaje('APLICAMOS TRANSFORMADA DE LAPLACE y subtituimos condiciones iniciales')
mensaje('SUBSTITUIMOS LA TRANSFORMADA DE LA ENTRADA')
edi=subs(edi,X(s), laplace(xi));
mensaje('DESPEJAMOS Y(s)')
%%% Para versiones superiores a 2016
mensaje('DESARROLLAMOS LAS FRACCIONES PARCIALES DE Y(s)')
%%%% Si se ejecuta en 2015 o menor comentar las 3 lineas anteriores
mensaje('Aplicamos transformada inversa, asi la solución es')
set(hFig, 'Position', [0 0 900 900])
axes1 = axes('Parent',hFig,'FontWeight','bold','FontSize',16);
fplot(xi,[0, t0],'b','LineWidth',2)
fplot(y,[0,t0],'r','LineWidth',2)
legend('Entrada x(t)','Salida y(t)','Location','Best')
xlabel('tiempo','FontWeight','bold','FontSize',16)
title('Entrada y Respuesta del sistema','FontWeight','bold','FontSize',16)
fplot(dy,[0,t0],'g','LineWidth',2)
title('Primera y segunda derivada de la salida','FontWeight','bold','FontSize',16)
fplot(ddy,[0,t0],'m','LineWidth',2)
legend('dy(t)/dt','d^2y(t)/d^2t','Location','Best')
xlabel('tiempo','FontWeight','bold','FontSize',16)